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A model representing the mechanical response of a rubberlike polymer is derived, using 
the same molecular concepts of entropy elasticity and a viscous parameter expressing 
interactions between polymer molecules employed by Rouse, Bueche, and Zimm in their 
molecular theories. Since the model developed here represents the mechanical response of 
a chain, rather than representing the chain itself, it can be modified more easily than these 
strict molecular theories to include effects due to entanglements between chains which 
modify the character of the viscosity-molecular weight relationship at a critical molecular 
weight. 

This modification is introduced, the results for both steady-state and transient response 
functions are calculated, and these results compared with experiment for the ''standard" 
polyisobutylene. The agreement indicates that the same ''entanglements," whatever their 
precise nature, are responsible for the proportionality of viscosity to M^-^ for a high molecular 
weight polymer, the steady-state elastic compUance, and the pseudo-e(|uilibriiini (•()m])hance 
at intermediate times or frequencies. 



1. Introduction 

Molecular theories of the viscoelastic behavior of 
rubberlike polymers have been developed in recent 
years by Rouse [1]\ Bueche [2], and Zimm [3]. In 
these theories the polymer molecule is represented by 
a model consisting of A^ segments connected by junc- 
tions which permit free rotation. The length of the 
segments is unspecified, but is assumed long enough 
to give a Gaussian distribution of end-to-end dis- 
tances, and hence account for an entropy type elas- 
ticity. The interactions between molecules are 
supposed to be purely viscous (for unvulcanized 
materials) and these viscous forces are lumped at 
the junctions as connections betw^ecn the typical 
molecule whose motions are analyzed and the aver- 
age medium, or in the case of the Rouse and Zimm 
theories the solvent. The Zimm theor}^ includes an 
approximate treatment of hydrodynamic interaction, 
that is the coupling of motions along a single chain 
through the solvent. 

Only a simple viscous interaction is included in the 
basic form of these theories. The entanglement 
effects which change the character of the viscosity- 
molecular weight relationship at some critical value 
of molecular weight have been included [4] by treat- 
ments which modify the calculated response func- 
tions rather than by including them directly in the 
original model. 

We are interested here in representing the me- 
chanical behavior of a high molecular weight, un- 
vulcanized, undiluted polymer in terms of a model 
whose elements have a clear relationship to molecular 
quantities. 'Ilie model we use is basically a phenom- 
enological one, since it represents the mechanical 
response of a polymer chain rather than representing 
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the chain itself. It does not seem feasible to include 
entanglement efTects in a model of the polymer mole- 
cule, but by dealing with a model representing 
mechanical response the required modification is 
relatively simple. 

The basic molecular theories have been sum- 
marized and compared in some detail by Ferry [5] 
in his recent book, so there is no need to review them 
here. The model presented in this paper was dis- 
cussed briefly in a previous paper [6] and its predic- 
tions for steady-state response given. We include 
here a comparison of the predictions of this model 
with those of the Rouse theory and the expressions 
for transient, as well as steady-state, behavior. 

2. Derivation of Model 

We present a derivation of the basic model, re- 
lating the elements appearing in a certain form of 
the response function to molecular concepts, and 
compare the predictions of this simple inodel wdth 
those of the Rouse theory. We next include a repre- 
sentation of the high frequency limiting shear com- 
pliance as a separate element, not connected with 
the entropy mechanisms dominant at lower fre- 
quencies. Mnally we modify the model to include 
the effects of entanglements characteristic of un- 
diluted polymers of high molecular weight. 

2.1. A Model Analogous to the Rouse Theory 

Assvune that the complex shear modulus of an 
unvulcanized polymer, (j*(^), where p is a complex 
frequency=6^+ico, can be represented by A^ terms of 
the form amP/(p4-^m)' This sum can be written as 
the ratio of two polynomials in p, both of order Nj 
pPn-i{p)IQn{p), where both P and Q contain 
nonzero constant terms. This ratio in turn can be 
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expressed as a continued fraction by a process of 
alternate division and inversion, giving for 



J*(^)-ai- 



1 



fiiV- 



«2 + 



1 



(2.1) 



where the a^ and /3's are functions of the X^'s. 

Since the q:'s in this expression represent a strain 
proportional to stress and the j3's a stress proportional 
to rate of strain, this response is represented exactly 
by the model of figure 1. 

This model represents the mechanical response of 
an unvulcanized polymer sample, and since the 
steady-state response is in fact one form of the 
general response function of a linear system it 
represents the response to a forcing function (either 
stress or strain) with any arbitrary time dependence. 
It is of the same form as the model used to represent 
a polymer molecule in the Rouse, Bueclie, and Zimm 
theories if we choose the constants as: 



l/ai = l/a2 = 



. =l/aN=NkT 



where N is the number of segments in the model 
representing the polymer molecule, n the number 
of molecules per cubic centimeter, and r\ the steady 
flow shear viscosity. Thus we should expect this 
model to give us the contribution of a single polymer 
molecule to the shear compliance of a bulk sample. 
We prefer to use constants that will give directly 
the response of a bulk sample, requiring the intro- 
duction of a factor n as shown in figure 1. 

The response represented by this model is the 
sum of responses of 1, 2, 3, . . . iV segments acting 
cooperatively. In this sense it corresponds directly 
to the normal mode calculation used in the develop- 
ment of the molecular theories. The assignment 
of constants follows the same molecular concepts 
used in formulating the model of the molecular 
theories, the springs representing the entrop3'-type 
elastic forces contributed by the segments of a 
chain, and the dashpots the viscous resistance 
opposing the motion of a segment. 



r^>-^ I 33-^ 



Figure 1. 



Simple model of mechanical response of polymer 
chain. 



2.2. Comparison With the Rouse Theory 

The steady-state response of the model of figure 1 
is given by : ^ 



J''{p) = l/(2NnkT) 

'^IprjnkTV' 



PV 



ANHkT 



)]' 



ctnh(A^a) (2.2) 



where 



sinh a 



=l\-^Ji- 



prj 



'N\nkTy^4.NHkT 

The Rouse theory gives : 
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(2.3) 



with 



_ 2N{N+2)nkT Jtt 



Both expressions depend on N, the arbitrary num- 
ber of segments into which the model of the molecule 
(or of the response) is divided. Rouse eliminates N 
from his results by approximations valid for frequen- 
cies less than a limit proportional to A^^^. Thus, to 
obtain expressions valid for all frequencies, we must 
take the limit as N goes to infinity. This is equiva- 
lent to assuming that segments of vanishing length 
can be treated as kinetic theory springs. 

Letting A^ go to infinity in the model used here 
eq (2.2) gives: 



G'^(j)), the reciprocal of J'^(p), is given by 



G*(p)^nkT ^-^tmh^^ 
or, equivalently, by: 



pri 
,kT 



(2.4) 



(2.5a) 



G^(2))=2nkT 22 -tV^ with X;;,=^ 

m=op-\-^m V 



(2.5b) 



The Rouse expression in this lunit is : 

6%p) =r-nkT JZ z^'^ X,-=-^^^$^ (2.6a) 



7 = 1 



P + \j 



677 



or 



G'(p) = (nkT/2)[y ctnh y-1]; 2/=^^^' (2.6b) 



Qprj 
Inkf 



These two functions are very similar, though not 
identical. Equation (2.5a), representing the response 
of our model, can be wi^itten as : 



2 This solution is obtained by the procedure described in L.A. Pipes, Applied 
Mathematics for Engineers and Physicists, pp. 265-267 (McGraw-Hill Book Co., 
NewTork, 1946). 
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G*(p) = (nkT/2)[x ctnh j'-x/sinh x] 



4 



x^-x -i-^ 



4pri 
nkT 



(2.5c) 



This is nearly the same as eq (2.6b), except that 
y = (3/2)^x. The second term in the brackets of 
eq (2.5c) becomes essentially unity for small values 
of l^^l, corresponding to that in eq (2.6b); for large 
values of \p\ the second term in both expressions is 
negligible compared to the first. 

2.3. Representation of the Limiting Compliance 

As noted in the last section, the process ol' going 
to the limit of infinite N in our model corresponds to 
taking segments of vanishing length, which results 
in expressions for J\cc) vanishing at high frequencies 
as w"^/^. We would not want to use a model which 
represented the limiting high frequency shear com- 
pliance as an entropy-type mechanism, since this 
would have an unrealistic temperature dependence. 
Therefore we go to the limit of vanishing segment 
length and represent this compliance by a separate 
element in the model. 

The model of figure 1 assumes that the junctions 
between segments are connected to the average 
medium by purely viscous eleinents, which show a 
zero compliance at infinite rrecjuency. In setting up 
this model we considered only the rc^sponse associated 
with rather long chain segments, and these viscous 
elements represent the resistance to motion experi- 
enced by the whole segment in its motion. In 
addition to this type of response, there also exists 
the type of compliance characteristic of polymeric 
glasses, which does not involve long segments of the 
chains but is connected with forces between indi- 
vidual atoms irrespective of whether or not they are 
a part of the same polymer molecule. This elastic 
mechanism will permit a small compliance, negligible 
at low frequencies in comparison to the entropy 
response, but becoming dominant as the frequency 
is increased. 

This type of response can be introduced into our 
model by connecting the viscous elements to ground 
through a very stiff spring, rather than directly. 
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Figure 2. Model including limiting high-frequency compliance. 



This modification is shown in figure 2, which we 
represent initially as a model with distributed rather 
than lumped parameters. The simplest assignment 
of constants, shown in the figure, is achieved by 
taking the total length of the model as unity and 
assigning the constants as shown to represent the 
response of the bulk sample. Ci is evaluated from 
the calculated response of the model at infinite 
frequency. In addition to the mechanical model, 
the electrical analog is also shown, using the system 
in which voltage corresponds to stress and current 
to rate of strain. 

2.4. Introduction of Entanglement Effects 

The model as formulated up to this point repre- 
sents the behavior of a collection of polymer chains 
in which the motion of any segment or group of 
segments will be opposed by a viscous force propor- 
tional to the length of that segment or group of 
segnumts. We know empirically [7] that this simple 
picture is not adequate to represent the viscous 
behavior of materials of high molecular weight, and 
that above some critical molecular weight, char- 
acteristic of each polymer, the steady flow viscosity 
becomes proportional to M^-'^, M the molecular 
weight. 

Bueche [8] has developed a theoretical picture in 
which this eft'ect is attributed to the efl'ects of 
entanglements between polymer chains, which be- 
come effective only when the chains reach some 
minimum length. Above this critical length he 
predicts that the viscosity should increase more 
rapidly than below it, going in the Ihnit of very high 
molecular weight to an expression proportional ta 
AP-K 

We adopt this theoretical picture, using however 
the empirical exponent of 3.4. We assume that the 
coordinated motion of any segment or group of 
segments with a length less than that corresponding 
to the critical molecular weight will l)e opposed by 
a viscous force proportional to the length of the 
segment or group of segments. When this length 
exceeds that corresponding to the critical molecular 
weight, Mc, we assume that the viscous resistance 
will still be proportional to length, but the propor- 
tionality constant will be much greater, representing 
the fact that the motion of such a section of the 
molecule will be opposed by its entanglements with 
other chains of the sample. 

The final model for the response, including 
entanglement effects, is shown in figure 3. Again 
this represents a model with distributed parameters. 
To simplify the final expressions, the length, I, of 
the model is defined as MjMc, with the length 
corresponding to the critical molecular weight taken 
as unity, r is defined as the viscosity of a polymer 
of molecular weight equal to the critical molecular 

weight; U is defined as ._ to make the viscosity 

represented by the model proportional to M^-^. c, 
corresponds to the equilibrium compliance of a 
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cross-linked polymer with molecular weight Mcl2 
between cross links, since this corresponds to the 
average molecular weight between entanglements. 

Although this is really a phenomenological model 
of response, the constants are simply visualized in 
terms of our molecular concepts. For motions 
connected with high-frequency response only short 
lengths of the polymer chain are involved, and hence 
only the resistive terms characteristic of small 
molecules are effective. For lower frequencies, 
where greater lengths of the chain are involved in 
cooperative motions, the elements representing en- 
tanglement effects become important. 

In analyzing the response of this model it will 
prove useful to employ various approximations valid 
in certain regions of time or frequency. For this 
purpose it is necessary to say something about the 
relative values of the constants appropriate for 
a model which represents the response of a rubber- 
like polymer. The constant C2, equal to the equilib- 
rium compliance of a vulcanized rubber with 
molecular weight Mc/2, will typically be of the 
order of 10~^ cmVdyne. r represents the viscosity 
of a polymer with molecular weight Mc, so should 
be within an order of magnitude of 10* poise. I is 
defined as MjMc] it must be greater than unity and 
in most cases where this model would be used it will 
be considerably greater. In the example considered 
at the end of this paper it is 60, yielding a value of 
i^=1.9X10*. {ciC2y''^ equals the limiting high fre- 
quency compliance which is of the order of 10~^^ for 
all polymers; thus Ci will be of the order of 10"-^*. 

3. Calculation of Steady-State Response 

The response of a polymer whose mechanical 
behavior is represented by the model of figure 2 
to any arbitrary (shear) stress or strain is found by 
imposing the given forcing function at the input 
and calculating the resulting response, again at the 
input. The stress required to maintain a constant 
rate of strain will, except for transient terms, be 
the integral over the line of the viscous term alone. 



This gives a viscosity equal to rl' 



.(M,(M.)' 



The steady-state response to sinusoidal stress or 
strain can be found almost as easily. This is the 
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FiGUKE 3. Model including effects of entanglements. 



problem of determining the input impedance of a 
two section electrical transmission line. From the 
differential equations relating the shear stress a to 
the shear strain 7: 



-2/(7 



(j.^ — ^7 



where 



'^=^ 



d7 



y=iooC2 



-ICCCiT 

rU 



l+iooCivL^ 



0<?<1 



we derive : 



--yzy. 



For steady-state sinusoidal stress and strain both 
force and response must be proportional to exp (icct) . 
Inserting this, plus the boundary conditions that 
a(l)=0 and that both stress and strain must be 
continuous at ^=1 we obtain: 



G*(ia)) = 



1 



_C7(0) 



J*(^co) 7(0) 

I^ tanh Ti + p tanh T^ 
Co 1 + p tanh Fi tanh F3 



(3.1) 



where 



-Vi 



io)C2r 



+ iooCiV 



F2 = 






ioiC2rU 






I? 



K=L{l-\) 



These expressions can be simplified immediately 
because of the fact that C2>>Ci and CirU<^<^ 1 for 
all cases of interest. Because of these relations, the 
denominators of Fi, T2, and r3 can be taken as unity 
for all frequencies up to those for which the hyper- 
bolic tangents become essentially unity. Thus, for 
real frequencies, eq (3.1) becomes: 
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G*{io^)=- 



■■yj'+i/ 



tanli a^L tanh Ka 
1+L tanh a tanh Ka 



W^+l' 



X 



(£+1) smh {K+l)a+{L-l) smh (ir-l)a: 
(L + 1) cosh {K+l)a-{L-l) cosh (i^-l)a 

a=^^'ia)C2r. (3.2) 



For our later extension of this expression into the 
complex plane to derive expressions for the transient 
response of the model, we shall have to go back to 
the full expression (3.1). 

Equation (3.2) can be coded for machine computa- 
tion in complex form, but it is still rather complicated 
for algebraic analysis. The relative values of tlie 
constants required to represent tlie behavior of a 
typical rubberhke polymer suggest several approxi- 
mations, valid in different frequency ranges. 

At very low frequencies the hyperbolic tangents 
can be expanded in power series, yielding for the 

low frequency limits ( oiC2rK^<^ 



its ( coC2riir^<^V 



G''{io^) = G\o^)+iG''{o,) 



co^C2r^ 



{LK'+ZUK'+'^LK+\)+ . 



+i[c^r{LK+l)+ .. .1 

o 

(foi-/'^»l). 
J*(icc)=J'ico)-iJ"(co) 

C2 LK'+SL^K'+'SLK+l , 



(3.3) 



-.[. 



{LK+ir 



o:r{LK+l) 



02 



(1+2)- 



ccrP 



(forP-^»l). 

(3.4) 



When CO becomes large enough so that tanh Ka, 
but not tanh a, is essentially unity, eq (3.2) reduces 
to: 

tanh a+Z 



a*(ic)=- 



C2 l+L tanh a 



=-ctnh (a+x); tanhx 
C2 



-r 



(3.5) 



shice ^ 1 + icoCir can still be taken as essentially unity. 
This can be separated into real and imaginary 
parts giving 

^, , . e sinh 2((9+x)+sin 2^ 
G M-=— 



G''{c.) = 



C2 cosh 2(6>+x) — cos 2(9 
e sinh2((9+x)-sin 2(9 



J'ic.)- 



C2 cosh 2 ((9 +X)— cos 2(9 

C2 sinh 2((9+x)+sin2(9 
'2^*cosh2((9+x) + cos2(9 



(3.6) 



T^w )—^ sinh 2((9+x)— sin 2(9 
^'^■^~29' cosh 2((9+x)+cos 2(9 



(3.7) 



where 



VcoC2r 



and 



25 25 



At still higher frequencies, both hyperbolic tan- 
gents can be taken as unity, but tlie denominator of 
Fi must be retained, giving us: 



0.=^ 



1 



VC1C2 L i 



ooCir 






(or co> 



25 



(3.8) 



The expressions for the imaginary parts of the 
modulus and compliance can be differentiated ex- 
plicitly and equated to zero to find the extrema of 
these curves. However, the resulting transcen- 
dental equations are too complicated to be solved 
other than numerically, so this procedure is not of 
much value. Certain approximate expressions for 
the extrema can be found, based on the fact that the 
functions are closely approximated by regions pro- 
portional to some power of frequency, with only 
narrow intermediate frequency regions. 

The imaginary part of eq (3.6) can be approxi- 
mated by tliree expressions, approximately valid in 
different frequency regions, and from the inter- 
section of these expressions approximate values for 
the maximum and minimum in G^^ obtained. In 
the region where both d and x must be considered, 
we obtain two approximations. 
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For^<X we find: G'' ' 



For 6>>x 



G''. 



C2X 

X 



The intersection of these two curves, at B'^=- x^ or 



1 



locates approximately the low-frequency 



maximum in G^\ 

When the frequency rises so that x can be neglected 

compared to 6, we predict G'^^-j and the intersec- 
tion of this with the second of the preceding two 
gives for a minimum: 



-p- or CO- 



4 
C2rL 



These expressions are fairly good approximations 
for the frequencies of the extrema of the calculated 
curves shown in figure 4. As might be expected, the 
corresponding predictions of the magnitude of G^^ 
at its extrema are not particularly good. 

By a similar procedure, the minimum and maxi- 
mum of J'^ are predicted. The frequency of the 
minimum in J^' is approximately 



that of the maximum 



or CO 



C2r V4i 



2/3 



2/3 



6'^ 



or 0)' 
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The high frequency maximum in G'^ can be ob- 
tained by the explicit differentiation of G^' in (3.8), 

yielding: .=^^. 
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Figure 4. — Corwparison of theory with experiment. 

Complex modulus, polvisobutylene. Lines, this theory; points, measurements 

[12]. 



4. Transient Response 

4.1. Theoretical Considerations 

The transient functions, G{t) the stress relaxation 
function, and J{t) the creep function, are given 
by [6]: 

G{t) = {ll2Ti) r°^'" {G'^ip)!!)) exp {pt)dp (4.1) 

J c-io3 

J«) = (l/2«) f'^'" {J*{p)lv) exp {j>t)dp (4.2) 

J C—ioo 

where (t*(^) and t/*(p) are the steady-state func- 
tions of section 3 continued into the complex fre- 
quency plane with p=s^ioi. It is convenient for 
subsequent use to write them in the form: 



G*{p)^llJ*{p) 

_T, (p+1) sinh (r3+r,) + (p-l) sinh {T,- 



■ro 



C2 (p+l) cosh (Tg 



ith 



-ri)-(p-i) cosh(r3-ri) 

(4.3) 



WiVr 



pCiV 



+pcir 



'^ti Vr 



pcivU 



-{-pcirU 



K 
L 



^3=^ r2, p=^' 






The functions G*(p)/p and J*{p)/p are both 
analytic throughout the complex plane except for 
certain singularities on the negative real axis, in- 
cluding one for J'*(p)/p at the origin, which is a 
second order pole. The others are first order poles 
except for two accumulation points at 

p' = -l/(cir) Sindp'' = -l/{cirU) 

which are essentially singular. In both cases the 
simple poles occur in two infinite sets, one between 
and p^\ the other between p^^ and p\ 

Both functions are single valued and both are 
even functions of Ti and r2. For this reason we are 
at liberty to choose either sign for the radicals, and 
we shall choose it to yield positive real parts for Ti 
and r2. 

No singularities occur to the left of ^'. The paths 
of integration can thus be deformed to a finite circle 
in the left half plane enclosing all poles except that 
at the origin for e/*(p)/p. This yields the following 
a priori bounds for G(t) and J{t): 

G{t)\<M, 

J(/)|<M2(l+0, M^ and M2 constants, 

for all ^>0. The term linear in t comes from the 
second order pole at the origin. 

Formal solutions cf (4.1) and (4.2) valid for all 
times may be written as: 
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J{t)- 



t 



r(l+LKyi+LK 



<?(0=2 9n exp (pj) (4.4) 

1 

-S jm exp (2)J) (4.5) 



C2 



where the poles pn and pm. are the roots of: 

(Pn + l) cosh (Ts,n+Ti,n) 

-(Pn-l) cosh (Ts,n-r^,n) = (4.4a) 

(p^+1) sinh (r3.^+ri. J 

+ (Pm-l) smh (r3,..-ri. J=0. (4.5a) 

From the rektions G^(0) = 6^*(oo) and J(0) = J*(oo) 
[9] we conclude: 



1 



?-''''=iw[^^+3|rqS)]-(^'^^^'' 



(4.6) 



(4.7) 



4.2. Expressions for the Long Time Range, 
^>10-W(/^+l)' 

For the very long time region tJie equations for the 
poles and corresponding expressions for the residues 
can be written as: 



(L+1) cos {K+\)r,-{L-\) cos (K-1)t„=0 

(4.8) 



9n 



2/C2 



l + (iiC/L)[l + (Z^-l)sin^rJ 



--{-PnC2ry 



and 6^(0 ~X) Qn exp {—Tlt/c2r). 



(i+1) sm (i^+l)M.+ (X-l) sin (K-1)m.=0, 



Mm=(-2>mC2r)^/^ 
(4.9) 



Jm 



2c,/m' 



l + (/^/L)[l + (L2-l)cos^M. 



and X; jm exp (^„0 ~S im exp i-1j.lt/c2r). 
1 1 

It is reiuHly seen that t„ and fi^ are bounded respec- 
tively by (n-l±^'^7r/(K+l) and (™±0 ^l(K+l). 
For the residues we have the bounds: 



(1 +KZ) - ' < C2j/„/2 < (1 +K/L) - ' 

These series are asymptotic in time, and in both 
cases the error involved in stopping at the Nt\i term 
will be less than 0.1 percent for 

t>c.fY{K+iyN-\ 

Our neglect of Ci in this region is justified since 
the value of PnCivU will be no greater than 
{N+l)V{c,lc2)L\K+l)-\ which is small for at 
least the first several hundred terms. In practice 
it is not feasible to compute more than 100 to 200 
poles and residues. Hence we turn to other methods 
for shorter times. 

4.3. Representation for the Short and Medium 
Time Ranges, /< 10 ' Co/'C/v+l)'- 

TJie series of the preceding section correspond 
to transforms of the Mittag-Leffier expansions of 
G*(p)/p and J*(p)/p. For shorter times we have to 
find expansions which are centered around infinity. 
Such expansions can be found in the form of expo- 
nential series in powers of exp ( — 2ri) and exp ( — 21^3). 
For convenience of notation we write them as Dirich- 
let series by introduchig tlie (k^finilions: 

Xi=exp (-2ri), ^3-exp (-2r3), f=(p-l)/(p+l). 
We find tlien: 



1_^ l±lL 



V 



J*(p)_ 



C2pl-^x, ^-x, 

1 Is 



1-f^l 



C2 1-fXi^ ""'l-fXi 






(4.10) 



(4.11) 



It can be shown that f, Xi, and x^ are less than 1 in 
absolute value for all points ofi' the negative real 
7;-axis. Hence there will be a region around the 

negative real axis outside of which ^ x^ will be 

less than one in absolute value. The following 
expansions hold: 



^^=I±j 1+2 S(m)' 



j^i+2f:(m)'} 

•{l-2..X3^±(^^.3y} (4.12) 
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1+2 2 (-r^: 

1 



X < l+2xia 









The series are asYinptotic in both Xi and x^,, i.e., if 
we terminate after n terms in x^, the remainder is 
0(2:3''+^). If we replace the second brackets of (4.12) 
and f4.13) merely by l+0(x3), it is not useful to go 
beyond the number of terms of the Xi expansion which 
results in a remainder term of ^{x^). 

We can estimate the error involved in truncating 
the Xi series by means of the following theorem [10]: 

Theorem: If L{x) is a function, continuous and 
positive on the real interval ()<x<Xq, with the 
property 

lim L(ux)/L(x) = l for all fixed 7/>0, 



then if y{F)^\f{p), with real F(t),is convergent for 
0<p<Po and/(^) has the asymptotic behavior 



andii F{f)>0, then 

C 



r;>o, 



«>o, 



r 



F(t)dt^ 



T{a+l) 



t-L(t) 



for ^->0. 



Furthermore, if F(t) is mono tonic in the neighbor- 
hood of ^=0, differentiation is allowed to yield the 
asymptotic behavior of F(t) itself. 

This theorem can be used to estimate the trun- 
cation error if we apply the inverse Laplace transform 
to a finite number of terms in eqs (4.12) and (4.13). 
We shall use it only to determine the range of 
usefulness of the leading terms, since these are in 
fact the only ones which can be transformed in 
closed form. Putting them on the left-hand side, 
we write: 



G{t)-G\t)^ 



^^^-M —2^Ux,y+o{xs)\ 






'^-' \ ^ fail 

I ClP 



} 



(4.14) 



(4.15) 



The expressions on the right under the i"~' signs 
satisfy the conditions of the tlieorem. After differ- 
entiation we get the following remainder terms: 



G{t)-CP{t) 

1 U-\ t 



■yJCiCi L"^ Ci?' 






1+1 



'•0^£f\ 



(4.16) 



J\t)-J{t) 



-^CxCi 




(4.17) 



G^it) and tF{f), the inverse transfonns of the lead- 
ing terms of eqs (3.12) and (3.13), are tabulated 
functions : 

G'{t)^-{c,c.^-"' exp (-fl2c,r)hit/2c,r) (4.18) 
J°(<) = (cic,)"^ e.xp i-t/2c,r) ^{l + t/c,r)h{t/2c,r) 

where Iq{x) and Ii(x) are niodifie<l Bessel functions. 
For very short times tliese may be expanded in 
power series, giving: 



(-1)'(2.)! 



J%t) 



G\t)=^[i+Y:, ,,,3 

VdCo L 1 ("0 

(-l)'+'(2.)! 



H'+?l^ 



\y(2p-i) V4c 






(4.20) 
(4.21) 



These expansions are useful in the range 0<^<10cir. 
For longer times we can employ as^nnptotic expan- 
sions of the Bessel functions which 3'ield: 



GV)^ 



jo(0~2^g[l-g:+...]. (4.22) 

We will be in error by less than 0.1 percent by 
taking G{t) = G^{t) for times less than C2r/182, and 
J(t)=J'^{t) for times less than C2^/49. We still need 
an expression to cover the range 

C2r/182<^< 10-V(i^+1)^ 

the upper limit being that established for the prac- 
tical use of the residue series. If we set Ci = in 
eqs (4.12) and (4.13) we obtain series which can be 
transposed term by term, giving: 
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where M< [K], that is we may not carry these series 
beyond [K] terms in view of our neglect of the x^ 
series {xz=xf for C]=0). 

The essential singularity at p^^ precludes any 
rigorous calculation of a remainder term associated 
with tlie neglect of Ci in this region of time. Widder's 
fornmla [11] gives a series of approximations of order 
A^ for G{t) in terms of ^*(2^) and its first A^ deriva- 
tives evaluated at p=N/t. For any finite approxi- 
mation, the terms in Ci vanish with increashig t, 
but the expression becomes exact only in the limit of 
N—^co, and no estimate is available for tlie error 
associated with tlie approximation involved in stop- 
ping with some finite N. 

A more convincing rationale for the neglect of Ci 
in this region can probably be based on physical 
considerations. Ci was introduced into the model 
to account for the glassy compliance, which nnist be 
added to the rubberlike compliance but is significant 
only at very short times or high frequencies. The 
steady-state functions change, in a fairly narrow 
frequency range, from the region where (j^ oc^V^ ^^^^ 
J' ccui~^f\ characteristic of the model with Ci = 0, to 
their limiting values of feci)"^^^ and {coCiY''^. At 
lower frequencies Ci does not appear in either G'^^iiS) 
or J*(ico). Since our expressions for G^{t) and J°(0 
are valid into the region of time where they are 
proportional to t~^^^ and ^^^^, corresponding to the 
transforms of p^'^ and p'^'"^, we should be quite 
justified in using the simpler model with Ci = at 
longer times. 

In practice we may need up to 250 terms in eqs 
(4.23) and (4.24), not a difficult task for a modern 
computer. We have purposely set the limits so that 
more terms are required here than in the residue 
series, since it is easier to compute the exponential 
and error functions than to determine the poles from 
the transcendental eqs (4.8) and (4.9). 

5. A Trial of the Theory 

To test the predictions of the theory, we utilize 
the expermiental results available on the NBS 
^ ^standard' ^ polyisobutylene.^ Although these re- 
sults are rather old, they still comprise the most 
complete set of data available for any single polymer. 
Direct measurements are available covering essen- 
tially the whole response curve for the steady state 
and stress relaxation functions, and creep measure- 
ments for times greater than and equal to that 
corresponding to the pseudo-equilibrium plateau. 



3 "Standard" here refers to a particular batch of polymer from which samples 
were distributed to permit comparison of different measurement techniques and 
results as described in reference 12. This comparison was concluded several years 
ago, and no further samples are available. 



The dynamic measurements used are those ob- 
tained by several investigators as summarized by 
Marvin [12]. The stress relaxation measurements 
are those reported by Tobolsky and (^atsiff [131, 
converted from Young^s to shear juodulus by the use 
of a factor )i throughout. Although this factor is 
not strictly correct near the glassy modulus, the 
error is not serious since the shear jnodulus apj)ar(Mit- 
ly never exceeds Ya of the bulk modulus. The creep 
data are unpublished results of Ijcaderman and 
Overberg using a concentric cylinder viscometei' [14]. 

The constants needed can be obtained from the 
empirical expressions for the viscosity of polyiso- 
butylene given by Fox et al. [7]. Values of /• and Cj 
were calculated from the critical molecular weight 
and the viscosity of polyisobutylene at this molecular 
weight. The measured viscosity and molecular 
weight of the ^ ^standard' ^ polyisobutylene are not 
consistent with the expressions given by Fox, 
probably because the milling required to obtain a 
large uniform batch produced a larger than normal 
low molecular weight fraction. Thus, for consistent 
results, we must either use a fictitious molecular 
weight corresponding to the measured viscosity or 
modify the constants of Fox's empirical relations. 

The hest rc^sults are obtained by the use of a 
fictitious juolecular weight of 1.03X10^, correspond- 
ing to the Jiieasured viscosity of tiie polyjner, rather 
than the value of 1.56X10^' from light scattering or 
even that of 1.35X10^ from intrinsic viscosity 
measurements. In a previous publication [0] the 
results of this same theory for the dynamic functions 
only were presented with / calculated from the 
relation r]=rP-\ the ezpoMent corresponding to 
Bueche's theoretical prediction [8]. Here we use 
'rj = rP-^^ corresponding to the empirical relation. 
The difference in tlie calculated results is small, but 
those obtained liere using the 3.4 power agree 
slightly better with the experiments than did the 
former' calculations. 

The constants used, calculated for a temper'ature 
of 25° C, and their source, are summarized as follows: 
c.2=MJ{2pRT) =3.77 X 10-^ M, from Fox et al [7] 
r==5.5X10*. Derived from the viscosity of polyiso- 
butylene with molecular weight of M^ [7] 
I =(r]/r)^^^-'^=Q0.l2, rj the measured viscosity at 
25 °C-10io-^l 



{'-f^']"" 



a .377X11)2 



if=i(/-l)-8.145X10'' 

Ci=(J'{o^)y/c^=l,S2XlO-'\ usuig the value J'(°o) 
from the dynamic measurements. 

The results are presented in figures 4, 5, 6, and 7, 
where the lines represent values calculated from the 
theory and the points the experimental measure- 
ments. The agreement is striking in view of the 
wide range of values covered. All the features of 
the experimental curves are reproduced qualitatively, 
with the changes in the character of the curves 
occurring at nearly tlie correct frequencies or times. 
The maximum deviation between the calculated 
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Figure 5. Comparison of theory with experiment. 
Complex compliance. 
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Figure 6. Comparison of theory with experiment. 
stress relaxation. 
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Figure 7. Comparison of theory with experiment. 
Creep. 

and experimental values is about 0.06 on a logarith- 
mic scale, correspondino" to about a 15 percent 
difference, except for a narrow region of frequency 
for J'' and G'\ where the difference reaches about 
0.3 logarithmically. The real parts of the dynamic 
functions and the transient functions agree more 
closely with the measured values. A result we did 
not expect is that the transient curves appear to 
be represented by the theory somewhat better 
than G^ and J\ 

The theor}^ is certainly not exact. In particular, 
the prediction of a proportionality to the square 
root of time or frequency (or its inverse) in the 



region between the rubbery and glassy plateaus is 
only approximated by the measured results. 

However, the theory yields quite satisfactory re- 
sults. They certainly agree with the measurements 
as closely as most calculations comparing different 
types of measurements through the use of approxi- 
mation methods now available. It is particularly 
interesting to note that the values of both the pseudo- 
equilibrium plateau and the equilibrium compliance 
agree with the measurements to well within the 
experimental uncertainties. These values as calcu- 
lated by the theory, it should be recalled, depend 
on constants obtained entirely from steady flow vis- 
cosity measurements. It is thus quite apparent that 
the same ^ ^entanglements '^ whatever their true na- 
ture, are responsible for the break in the viscosity- 
molecular weight relation and thepseudo-equilibrium 
rubbery plateau, as assumed in the theory. 

The necessity, in these calculations for the '^stand- 
ard'' polyisobutylene, of using a fictitious molecular 
weight precludes any chance of using these results 
to check on the effect of molecular weight distri- 
bution. This question must remain in abeyance 
pending calculations for other s\^stems, preferably 
fractionated polymers. 
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